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Duality Relations and Exotic Orders in Electronic Ladder Systems 
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We discuss duality relations in correlated electronic ladder systems to clarify mutual relations 
between various conventional and unconventional phases. For the generalized two-leg Hubbard 
ladder, we find two exact duality relations, and also one asymptotic relation which holds in 
the low-energy regime. These duality relations show that unconventional (exotic) density-wave 
orders such as staggered flux or circulating spin-current are directly mapped to conventional 
density-wave orders, which establishes the appearance of various exotic states with time-reversal 
and/or spin symmetry breaking. We also study duality relations in the SO (5) symmetry that 
was proposed to unify antiferromagnetism and d-wave superconductivity. We show that the 
same SO(5) symmetry also unifies circulating spin current order and s-wave superconductivity. 

KEYWORDS: Electronic ladder, Hubbard model, duality relation, exotic order, S0(5) symmetry, staggered 
flux, spin circulating current. 



1. Introduction 

Unconventional density-wave orders such as <i-density 
wave (dDW) (which is equivalently called as staggered 
flux or orbital antiferromagnets) and (i-spin-density wave 
(dSDW) (which is circulating spin current) were first pro- 
posed in the context of excitonic insulators 1 and later 
discussed in high-T c superconductors, 2-4 but the appear- 
ance of these orders was not established at that time. Re- 
cently, several experimental results have led to a resur- 
gence of interest in the possibility of these exotic orders. 
A dDW state 5 was discussed to appear in the under- 
doped region of high-T c superconductors, where a pseu- 
dogap was observed, 6 and also in the low-temperature 
phase of the quasi-two-dimensional organic conductor, 9 
a-(BEDT-TTF) 2 KHg(SCN) 4 . Both a dSDW state 7 and 
a dDW state 8 were proposed as origins of hidden order 
in the hcavy-fcrmion compounds URu2Si2 and UPt3. 

The appearance of unconventional (exotic) orders has 
been tested in microscopic models of correlated elec- 
trons. In particular, the generalized Hubbard model on 
the two-leg ladder has been attracting attention as a 
minimal model for showing the exotic orders. 10-15 As is 
well known, the standard two-leg Hubbard and t-J lad- 
ders without any inter-site interaction do not show any 
order. 11 ' 14 ' 16 Both strong-coupling and weak-coupling 
analyses however found that the generalized Hubbard 
ladder model with inter-site interactions exhibits var- 
ious phases at half-filling. 10 ' 12 ' 13 ' 15 There are at least 
eight phases, 13 i.e., charge-density-wave (CDW), dDW, 
p-density wave (pDW), /-density wave (/DW), $• 5- 
Mott, and )-Mott phases. For less than half-filling, 
large-scale numerical calculations also reported CDW 17 
and <£DW 14 phases. Bosonaization and renormalization- 
group analysis revealed the appearance of quasi-long- 
range order of density waves. 15 ' 18 Recently we found two 
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exact duality relations between these various phases. 19 
The relations show that unconventional density- wave or- 
ders such as staggered flux or circulating spin current are 
dual to conventional density-wave orders and there are 
direct one-to-one mappings between dual phases in the 
generalized Hubbard ladder systems. 

The electronic ladder system also serves as a play- 
ground for theories of high-T c superconductivity. The 
SO (5) theory, which was proposed to unify antiferro- 
magnetism (AFM) and d-w&ve superconductivity (eJSC) 
in terms of the SO (5) symmetry, 20-22 was also applied 
to the Hubbard ladder. Scalapino et al. 23 and later 
many groups 10-12, 24 studied the Hubbard ladder with 
the SO (5) symmetry to capture some of the basic low- 
energy physics of the high-T c cuprates. There, one grand 
order parameter field 



(1) 



behaves as a five-component vector, where Odsc denotes 
the pairing operator of dSC with ReO = \{0^ + O), 
ImO = — O), and the three elements Af a with 

a = x,y,z are Cartesian components of the staggered 
magnetization. 

In this paper we develop further the duality relations 
given in ref. 19 and also study the duality structure in the 
SO(5) symmetric Hubbard ladder. First, we discuss two 
exact duality relations in the generalized Hubbard lad- 
der. One duality transformation relates conventional den- 
sity waves to unconventional density (or current) waves 
and the other relates charge-density degrees of freedom 
to spin-current ones. These transformations also show 
duality relations between s- and d-wave superconduc- 
tivity (sSC and dSC). Furthermore, we show another 
asymptotic duality relation, which appears only in the 
low-energy effective theory at half-filling. These duality 
relations among various density-wave phases are sum- 
marized in Fig. 1. The transformations give one-to-one 
parameter mappings between dual phases. If one finds an 
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Fig. 1. Two exact duality relations (I and II) between four phases 
with (a) CDW, SDW, dDW, and dSDW orders, and (b) pDW, 
/DW, pSDW, and /SDW orders. There is also another duality 
relation (III), which appears only in the low-energy region at 
half-filling. 



ordered phase in a certain parameter space, one can nat- 
urally conclude the appearance of dual ordered phases 
in dual parameter spaces. These relations help us to find 
various new exotic phases with spin and/or time-reversal 
symmetry breaking. 

Next, we apply the duality relation to the SO(5) sym- 
metry. We find that the SO (5) symmetry also unifies 
dSDW (circulating spin current) order and s-wave SC. 
Also, the SO (5) symmetric Hubbard ladder has a dual- 
ity structure. In one parameter region, the original AFM- 
g?SC SO (5) vectors represent dominant correlations, but 
in the other dual region the new dSDW and sSC vectors 
are dominant showing a crossover from dSDW to sSC 
upon doping. 

This paper is structured as follows. The definitions 
of the Hamiltonian and operators are given in Sec. 2. 
The two exact duality relations I and II are described in 
Sec. 3. The duality structure in the generalized Hubbard 
ladder model is shown in Sec. 4. Section 5 contains a 
discussion on the duality in the SO (5) symmetry Dual- 
ity relations are reconsidered by using the Bosonization 
framework in Sec. 6. An application of the duality rela- 
tions to results in the strong-coupling limit is discussed 
in Sec. 7. Section 8 contains discussions. 

2. Definitions 

2. 1 Hamiltonian 

We consider the generalized two-leg Hubbard ladder, 
which contains the on-site repulsion U , the intra-rung 
repulsion Vj_, the intra-rung spin-exchange J± (Jf), and 
the intra-rung pair hopping i pa ; r . The Hamiltonian is 
given by 



H = Hq + H It 



(2) 



with 



Ho 



t±Cj thtr C j>2 ,c r } + H.C.], 



Cj + 1,1,<t + cj :2l( T C J + l,2,<x) 



(3a) 



fining = 53 ^ 53 n iM n M + ^ 53 
3 L I cr.cr' 

+ *pair(ct 1)T ct 1)i C,-,2,lCj >2 ,t + H.C.) 



(3b) 



Here c ].i.cr ( c 3-i^) (I = 1, 2 and a =T,I) denotes 
an electron creation (annihilation) operator on the lih 
site of the jth rung, rij t i t(T — ct l „Cjj,, a an d Sj.i — 

\ Eo- a' c ] i <j CT < J °' c iJ-,' J ' w ith the Pauli matrices a a (a = 
x,y,z). In this paper we consider both spin isotropic 
(J± = Jf) and anisotropic cases. 

2.2 Order operators and typical states 

Here, we list the definition of order operators used 
in the present work. We consider a charge-density-wave 
(CDW) operator (or s-density-wave operator) 



a 



CDW 



(i) = l^2(nj,i,* 



(4) 



a ci-density-wave (dDW) operator 

OdDwti) = \ E(c],i,<t c j,2,<t - H.c), (5) 



a p-density-wave (pDW) operator (or staggered dimer 
operator) 



O 



pDW 



(■?') = I H[( C I + 1,1,<t C JM^ - C ]+l,2 l0 - c J,2 )0 -) + H.C.], 



(6) 

and an /-density-wave (/DW) operator (or diagonal cur- 
rent operator) 

OfDw(j) = J Y^4+l,l,<r C ^" ~ C ]+1,2,o- c jM,o-) - H .C.]. 

(7) 

In the similar way, we consider operators in the spin 
sector. Inserting the Pauli matrix a^ a into the right- 
hand sides of eqs. (4)-(7), one defines a spin-density-wave 
(SDW) operator 



OSDWC?) = - 53 a *<r( n 3,l,<r ~ n 3,2,<r)> 



(8) 



a d-spin-density-wave (dSDW) operator 

OdSDwti) = r 53 a ™( C j,hor C 3,2,<' ~ H - C - 



(9) 



a p-spin-density wave (pSDW) operator O p sdw(j) 



4+i,2,c7 C i,2,<r) + H.c], and an 



/-spin-density wave (/SDW) operator 0/sdw(j) = 

lE ff C[( c |+l,i,^-V ~ c ]+i,2,<t c j,i^) + H.c], respec- 
tively. Note that the CDW and pDW (SDW and pSDW) 
orders are kinds of density waves of charges (spins) while 
the dDW and /DW (dSDW and /SDW) orders finite 
local charge (spin) currents. Accurately speaking, e?DW 
and /DW operators do not correspond to exact current 
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operators if there are pair-hopping terms (i pa ir 7^ 0) in 
the Hamiltonian, but they can detect time-reversal sym- 
metry breaking, which is closely related to currents. Or- 
der parameters on the ladder are given by 

A (q) = L- 1 Y / OA(j)exp(iq 3 ), 

3 

where L is the number of rungs. 

Moreover, we consider the d- and s-wave pairing oper- 
ators, 

Cdsc(j) = -^=(cj- i iT Cj- 2,1 - cj. ld c jt 2 ,t), (10) 

O s ScU) = ^( c j,l,T c J,l,J. + c j,2,TCj,2,i), (11) 

which characterize the d- and s-wave spin-singlet super- 
conductivity (dSC and sSC), respectively. Using these 
operators, the representatives of the d- and s-wave Mott- 
insulating (d- and s-Mott) states at half-filling are given 
by 

|d-Mott) = n°LcC?')|0), k-Mott) = n°lsc0')|0), 



respectively, where |0) is the vacuum of the electron op- 
erators. We also consider the spin-triplet and spin-singlet 
pairing operators with odd parity given by 

Ot,oU) = ^(cjm.tcjai + cj,i,i c i.2,r). ( 12 ) 

°s,o{j) = ^(Cj.ljCj, 1,1 - Cj,2,TCj,2,i)- (13) 

3. Exact duality transformations I and II 

Here we describe the duality transformations on elec- 
tron operators 19 so that this paper is self-contained. 
They are given by gauge transformations on the bonding 
and antibonding operators dj,±,cr = { c j,i,a ± c j,2,a) / V%- 
Two duality transformations are presented in the follow- 
ing. 

3.1 Duality relation I: density and current 

Consider a gauge transformation of antibonding oper- 
ators given by the unitary operator 



U z (0) = Y[eM-i0d],-,^-^- 



(14) 



In the case of — ir/2, this operator gives the dual- 
ity transformation I, dj,± l<7 = Ui(ir /2)dj t ± ta Ui(ir/2)~ 1 , 
which yields 

d j,+,a = dj + a: dj—jj = idj—jj (15) 

for a =t, {. In terms of the electron operators Cj t i t<7 , this 
transformation is written as 



(16) 



Applying the transformation (15), we can obtain dual 
representations of operators straightforwardly as follows: 
The density-wave operators are transformed as 



Ccdw = —OdDW, 
CsDW = — CdSDW, 
CpDW = — C/DW) 



OdDW = CCDW, 
CdSDW = CsDW, 
C/DW = CpDW, 



(17) 



CpSDW — ~ O/SDWj C/SDW = CpSDW- 

Thus, current-order operators, such as dDW, dSDW, 
,/DW, and /SDW operators, turn out to be dual to 
density-wave operators, CDW, SDW, pDW, and pSDW 
operators, respectively. It should be noted that the uni- 
tary operator Uj{0) gives a continuous transformation 
between dual operators. For example, if one considers 
CDW and dDW operators, they are transformed as 

Ui(9)OcdwU i (9)- 1 = OcDwcosfl - O dD w sin 0. (18) 

The s- and d-wave pairing operators are also related 
by the transformation (15) as 



CsSC — CdSC- 



(19) 



Thus the d-wave SC phase is dual to the s-wave SC phase. 
By definitions, d-Mott and s-Mott phases are also dual 
to each other, 

|d^Mott) = | s-Mott), (20) 

where we omitted a constant factor. On the other hand, 
the parity-odd pairing operators (12) and (13) are invari- 
ant under the transformation. 

3.2 Duality relation II: density and spin current 

Consider a gauge transformation given by the unitary 
operator 

Un(O) ;«<. .,/,..■ ■ d] d,. tl )}. (21) 

3 

This unitary with = ir/2 gives the duality transforma- 
tion II, dj,± 1<7 = Uii(ir/2)dj_± M Uii(ir/2)~ 1 , which yields 



d 3-A = d j,-,h 



(22) 



dj,+,i - dj,+,i, 

d 3~,l = id 3-A- 

In terms of the electron operators Cj t i t(T , this transforma- 
tion is written as 



Cj,l,t 


= (e^M.T + 


.e 3 ™/ 4 c,, 2 , T )/V2, 






= (^i,t " 


f e"/ 4 c,- 2 , T )/V2, 






= (e W S,U + 


■e-^c jM )/V2, 






= (e-' i/4 C i, U 


+ e^% 2d )/V2. 


(23) 



It is easily shown that the transformation (22) gives 
duality relations between density-waves operators as fol- 
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lows 



OCDW 


= OdSDW, 


OdSDW 


= -0CDW 




= -0SDW, 


£*SDW 


= OdDW, 


Opuw 


= C/SDW, 


O/SDW 


= -CpDW 


C/DW 


= -CpSDW, 


CpSDW 


= C/DW- 



(24) 

Density waves of charges are transformed into spin cur- 
rents while density waves of spins into charge currents. 
This transformation thus exchanges density and current 
as well as spin and charge degrees of freedom. Similarly 
to Ui(6), the unitary operator Uu(6) gives a continuous 
transformation between dual operators such as 

Un^OcDwUnie)- 1 = CDW cos9 + O dSDW sm9. 

(25) 

We note that under the transformation (22) the d- 
wave and s-wave pairing operators, Odsc and O s sc, are 
invariant except for phase factors, and hence the d- and 
s-Mott states are also invariant. On the other hand, the 
triplet parity-odd paring operator £> t)0 is converted to 
the singlet parity-odd one Si0 , 



Ot,o — Os,o- 



(26) 



4. Duality relations in the electronic ladder 

Here, we apply the duality transformations to the 
Hamiltonian (2). It can be shown that these transforma- 
tions map the model onto the same Hubbard ladder with 
different coupling parameters. It is convenient to rewrite 
the Hamiltonian in terms of bonding and antibonding 
operators in the forms 

H ° = - E [*H E (4a,^+i,a, CT + H.c.) 

+ ^(4+,<A+,a - 4- l( A-, CT )] , (27) 

tfrung = [ A ( d ],+A d i <->l d j,+,l d i ' -'J- + H ' C -) 
3 

+ B(4 +(T d J -,_ lT 4_ ii d J - 1+ ,i + H.c.) + C^nf^nf^ 

+^E»SS,t»SSa+^E»53,t«5f-M]. ( 28 ) 

A=± A=± 

where n^jj a = dj A a dj^ a , and A = + (— ) represents the 
bonding (antibonding) orbital. The coupling constants in 
cq. (28) are given by 

A = (U - V± + Vir)/2 + (2Jl + JD/8, 

B = (U-V±- tpair)/ 2 - (2J± - J±)/S, 

C = V^ + Ji/4, 

D = ((7 + V± + tpair)/2 - (2J.L + Jl)/8, 
E=(U + V±- tpair)/ 2 + (2J± - JD/8- 

Under the transformations (15) and (22), the total 
charge density, total magnetization, and kinetic energy 



term Hq are invariant. Hence the chemical potential /x, 
magnetic field, and the parameters tu and t± are un- 
changed through the mapping. On the other hand, the 
intra-rung coupling terms in H rung are mixed up by the 
transformations. The parameter mappings are given in 
the following. 

4-1 Duality relation I 

It is easy to see in eq. (28) that the transformation (15) 
changes only the sign of the A-term, but keeps the rest of 
terms invariant. Hence, it gives a one-to-one parameter 
mapping 



(A, B, C, D, E) -» {-A, B, C, D, E). 



(29) 



This leads to an exact duality relation in the parameter 
space of the system. The model with a parameter A is 
dual to the model with —A and self-dual in the space 
A = 0, i.e., 

U-V±+ tpair + (2 J± + Ji)/4 - 0. (30) 
The coupling parameters are mapped as follows: 

U=(U + V±- tpair)/2 - (2J ± + JD/8, 

V± = (U + W ± + t pair )/4 + (2J ± + JD/16, 

tpair = {-U + V ± + ipair)/ 2 - (2J_L + J Z ± )/8, (31) 
J± = -U+V±- tpair + (2J± - Jl)/4, 
Jl = -U+V±- ipair - {2J± - 3JD/4. 

Note that the spin isotropy (Jj_ = JJ) is conserved 
through this parameter mapping, i.e., J± — JJ. 

From the mapping (31) and the duality relation (17), 
one can conclude that if a density- wave order, e.g., CDW 
or SDW order, appears in a certain parameter region, 
a dual current order, i.e., c?DW or dSDW order, re- 
spectively, exists in a corresponding dual parameter re- 
gion. Because of this duality relation, all phase bound- 
aries must be symmetric with respect to the self-dual 
space. Indeed, the transition line between the CDW and 
JDW phases derived for half-filling and in the weak- and 
strong-coupling limits 10 ' 12, 13 coincides with the self-dual 
line (30). We stress that our exact result holds in gen- 
eral cases, regardless of the coupling strength, filling, and 
system size. 

4-2 Duality relation II 

From eq. (28), one can see that the transformation (22) 
changes only the sign of the B-term and hence gives a 
one-to-one parameter mapping 



(A, B, C, D, E) -» (A, -B, C, D, E). 



(32) 



Thus the present model has another duality: The model 
with a parameter B is dual to the model with — B and 
self-dual in the space B — 0, i.e., 



U-V±- t pair - (2J ± - JD/4 = 0. 



(33) 
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The transformed coupling parameters are given by 
U = (U + V± + tpair)/2 + (2Jj_ - J x )/8, 
Vx = (U + W ± - t pair )/4 - (2Jl - JD/16, 

tpair = (U - V± + tpair)/2 - (2Jj_ - J|)/8, (34) 
J X = C/ - F ± - i pair + (2J1 + JD/4, 

J* = -C/ + y ± + t pair + (2J± + 3JD/4. 

This duality relation leads to the conclusion that, if 
CDW or e?DW order, for example, appears in a certain 
parameter region, spin current (dSDW) or SDW order 
exists in a dual region, respectively. Note that even if we 
start from a spin isotropic model (Jj_ = JJ), the dual 
model is spin anisotropic and hence the spin symmetry 
breaking associated with SDW or dSDW order can occur 
in the dual model. The whole phase diagram must be 
symmetric with respect to the self-dual space (33) and 
the direct phase transitions between dual phases, if ever, 
locate exactly on the self-dual space. 

4-3 Duality relation between spin and charge 

A combination of the duality transformations (15) 
and (22) leads to the spin-charge duality transformation 
given by 

d J,+,l = d j,+,h = — *^J>+>J.> (35) 

d j ,-,T = 4? ,-,T> d j,-,l = id j,-,l> 

which directly exchanges charge and spin de- 
grees of freedom. For example, general density- 
wave operators J2kv /i(k)4 (k)c CT (k + Q) are 
transformed to spin-density-wave operators 
EkE<7 )£ r'=r,| /A(k)4(kK CT ,cv(k + Q). 

The transformation (35) gives another duality relation 
in the Hamiltonian, which is given by the combination 
of eqs. (31) and (34). The model is self-dual in the space 
satisfying both U - V± + JJ/4 = and t pail + Jj_/2 = 0, 
which is the intersection of the self-dual spaces (30) for 
duality I and (33) for duality II. 

4-4 Quantum phase transition at self- dual points 

Here we briefly discuss the nature of quantum phase 
transitions at self-dual points. In the self-dual spaces 
(30) and (33), the model Hamiltonian is invariant under 
the continuous rotations given by the unitary operators 




Fig. 2. Orders of dual phases are destabilized at the self-dual 
points, i.e., A = for the duality I and B = for the dual- 
ity II, because of the hidden U(l) symmetries. Hence the direct 
phase transitions between dual phases are of second order. 



Uj(0) and Uu(6), respectively. [This can be easily seen 
in eq. (28) by setting A = or B — 0.] Hence the self- 
dual models have extra hidden U(l) symmetries. Because 
of the U(l) symmetries in the self-dual space, a rigorous 
theorem 25 concludes that the dual orders that are contin- 
uously transformed to each other disappear on the self- 
dual models in one dimension. Hence, the direct phase 
transition between the dual phases, if it exists, must be of 
second order. See Fig. 2. An exception can appear only if 
susceptibility of the generator for this rotation is diverg- 
ing at the transition (see ref. 25); In this case, the orders 
can survive even at the transition point in the self-dual 
space, but there are gapless excitations associated with 
the continuous U(l) symmetry breakdown. We note that 
the discussion above does not exclude the presence of 
gapful phases in the self-dual models. Self-dual phases, 
which can be either gapful or gapless, may appear in a 
finite parameter region including the self-dual space. 

In the quantum phase transitions at self-dual points, 
the A- and .B-terms in the interaction (28) serve as 
symmetry-breaking perturbations for the U(l) symme- 
tries associated with Uj(6) and Uu(6), respectively. Ac- 
tually, the A- and £?-terms can be expressed with the 
difference of dual operators in the forms 

A i d ),+A d o--A d ],+.i d 3,-A + H - c -) 
= A[(O cdw ) 2 - (O dDW ) 2 } - A[(O dSDW ) 2 - (Osdw) 2 ] 
= A[\O sSC \ 2 ~ \O dSC \ 2 }, (36) 

B(d] +/t dj _ jT dt _ :l d ji+il + H.c.) 
= ^[(Ocdw) 2 - (O dSDW ) 2 } - B[{O c mw) 2 - (Osdw) 2 ] 

= B[\O s A 2 -\OtA 2 l (37) 

where \0\ 2 = (RcO) 2 + (ImO) 2 . When a symmetry- 
breaking perturbation is relevant, it induces a order and 
a gap determining the criticality of the phase transition 
(Fig. 2). For example, the transition between CDW and 
dDW phases was studied using the Bosonization method 
in the weak-coupling limit and shown to be characterized 
by c = 1 Gaussian criticality. 10, 13 

5. The SO(5) symmetric ladder 

The SO(5) symmetric ladder presented by Scalapino 
et al. 23 has only intra-rung couplings and belongs to the 
Hamiltonian (2). They showed that the half-filled Hub- 
bard ladder (ju = |J7 + V) has the SO (5) symmetry if 
the parameters satisfy 

Jx = Jl =4(U + V±). (38) 

No condition is imposed on the values of t±, ty , and i pa i r , 
since the hopping terms in Hq and the pair-hopping term 
(ipair) are SO(5) symmetric. 23,24 In this section, we show 
a duality relation in the SO(5) symmetric ladder. 

Scalapino et al. 23 showed that the five-dimensional 
SO(5) superspin vector n a (a — 1,...,5) is related to 
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the AFM and dSC operators by 



n 1 = V2ReOdsc, 
n 5 = V2ImO dsc , 



n 



(2,3,4) 



(x,V,z) 



?{x,y,z) 



(39) 



where we have omitted rung indices j. In terms of 
the rung charge Q = |^ ff ( n i<i + n 2<r - 1), the rung 
spin S a = Sf + , and the n a operators — 

~|Z)<7,<7' c L( (J a cr y)<T<T' c L'' tnc ten-dimensional SO(5) 
symmetry generators i afc (a, 6 = 1, . . . , 5) are expressed 
with 



Since the set of SO(5) generators are invariant under 
the duality transformation I, the SO(5) symmetric Hub- 
bard ladder is mapped onto the SO (5) symmetric model 
satisfying the condition (38). One can easily check that 
the condition (38) is conserved under the parameter map- 
ping (31). However, the parameter point in (Jj_, U, t pa i r ) 
for the SO(5) symmetric model is converted to the dual 
point with 














+ 







7rt 


+ 7Ty 


-s z 





f 

"z 


+ TTz 


sy 


-s x 




Q 









(4-«z) o 

(40) 

where the elements are antisymmetric. The SO (5) scalar 
operator p is expressed by the charge density operator 
with 

P = Ocdw + 1. (41) 

Applying the duality transformation I, we find that 
the transformed superspin vector h a is related to the 
spin current (dSDW) and sSC operators by 

h 1 = VmeOssc, n 5 = V2lmO sS c, 



fi (2,3,4) 



j 



(x,y,z) 



(42) 



where the three elements jf (a — x,y, z) are Cartesian 
components of the spin current (dSDW) along each rung 
and j z = OdSDW- Moreover, using the relations S = S, 
7r = iir, and Q — Q, one can find that the transformed 
symmetry generators L ab are given by just a permutation 
of original generators L ab , 

£(2,3,4)1 = £5(2,3,4) £5(2,3,4) = _£(2,3,4)1 (43) 

and the rest of L ab (a > b) arc unchanged. Thus, the 
set of SO (5) generators are invariant under the duality 
transformation I, except for the permutation of elements. 
From these results, we can conclude that, besides the 
superspin vector (39), the same SO(5) symmetry with 
the generators (40) unifies the dSDW and sSC in another 
grand order parameter field 

(V2lmO sSC ,f s ,j y s ,fs, -V2ReO sSC ). (44) 

The transformed SO (5) scalar operator p is given by the 
dDW operator with 



p = -CdDW + 1- 



(45) 



One may be surprised that the dDW operator is an SO (5) 
scalar, but one can easily express the dDW (and CDW) 
operator in terms of SO(5) spinors in an SO(5) symmet- 
ric form. 



U — —U — V± — ip a i r /2, J± — 2V± — tpair, 



t 



pair 



-2U -V±+ ipair/2, 



(46) 



and V± is given by V± = jj_/4 — U. 

To demonstrate the duality relation in the SO (5) sym- 
'metric model, we write down the Hamiltonian in terms 
of dual and self-dual operators. The intra-rung part (3b) 
can be cast into the form 



H r , 



+ 



3,1," 3 



J± 

4 



t 



pair 

2~ 



£(^f) 2 



a<b 



tpair (pj - 1) - 2tpair(/5j - 1)" 



(47) 



where no = \U + V. The Casimir operator C = 
^2 a< b(L ab ) 2 is invariant under the duality transforma- 
tion, and the kinetic energy term as well. It is also in- 
structive to rewrite the Hamiltonian in terms of the dual 
superspin vectors in the form 



"i £ ' 

3,1," 



H Iung — no n 



= £ 



20 



AU 7t r _ 



10 



£(^f) 2 



(48) 



a<b 



- 1 (J ± + AU - 2t p air) EK) 2 + £(*$) 

a a 

where we have used the Fierz identity 23 

a<b 



(49) 



The duality relation in the Hamiltonian is transparent in 
the forms (47) and (48). The SO (5) symmetric Hamilto- 
nian is self dual if two coefficients of dual operators are 
equal, i.e., 



J± + AU + 2t p air = 0. 



(50) 



The self-dual SO (5) model has a superspin SO (5) x du- 
ality U(l) symmetry, in total. The two-dimensional self- 
dual plane divides the three-dimensional SO (5) symmet- 
ric parameter space into two subspaces. In one subspace, 
the AFM and dSC correlations are dominant, showing a 
crossover upon doping as was discussed in refs. 20-23. In 
the other subspace, the <iSDW and sSC correlations come 
out to be relevant orders and hence symmetry breaking 
perturbations enhance dSDW or sSC order in the same 
way as the case of the AFM and dSC orders. 
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For tp a i r = and U < 0, this self-dual plane (50) 
coincides with a phase boundary line found in the pre- 
vious studies. 10 ~ 12 < 23 > 24 For example, the transition line 
between the CDW and dDW phases given in ref. 11 is 
identical with the self-dual line. For £ pa j r = and U > 0, 
on the other hand, the self-dual space locates in the C2S2 
critical phase. This implies that the critical phase is self- 
dual and a crossover between dual states characterized 
by different dominant correlations occurs on this self- 
dual line. 

6. Duality relations in weak coupling 

Here we reconsider the duality relations in the weak- 
coupling limit. In this limit, the duality relation I cor- 
responds to an already known relation in Bosonization 
studies. 10, 12, 13 We also find another new duality relation 
which appears only in the low-energy effective theory. 

6.1 Bosonization framework 

We apply the Abelian bosonization method. 26,27 In a 
continuum limit, we expand the electron-field operators 
as 

^ Xr7 (x) = e ik ^ x ^ RXa (x) + e- ik ^ x i> LXa {x) (51) 

for A = ± and a =T,|, where ipR\ a (ipLXa) represents 
the chiral fields for right- (left-) moving electrons in the 
bonding (A = +) and antibonding (A = — ) bands and 
kp\ are the Fermi wave vectors. As in the usual way, we 
introduce boson fields (p p xa for the chiral fields as 

^p\a(x) = J^Z—. exp[iptp pX(7 (x)}, (52) 

ior p = R/L = +/— , where n Xa denote the Klein factors 
and ao the lattice constant. Then, we define a new set of 
boson fields: 

<f>vr = (f>„ r + 4\, r , 6 Vr = (f)^ r - (j)^ r , (53) 

with 

€r = \{¥>p+1 + fp+l + r (<Pp-l + <P P -l)}, (54a) 

C = \Wp+] - Vp+l + r{<p P -] - <P P -l)} (54b) 

for v — c,s and r = ±. The <fi and 6 fields (53) 
satisfy the commutation relations [(j) vr (x),9 l/ i r i(x')] — 
—iir<d(x' — x)8 vl/ >5 rr t [0(x) is the Heaviside step func- 
tion] and [4> vr {x),(j> v ' r i{x')] = [0 ur (x),6 v > r >(x')] = 0. 

Using these boson fields, one finds that the duality 
transformation I is expressed as the translation of the 
6 C r fields by rir/2 while the duality transformation II is 
the translation of the 6 C + and S - fields by n/2. It was 
already realized that the former transformation gives a 
duality relation in the Bosonization method. 10,12 ' 13 

6.2 Duality relation III 

Here we restrict our discussion to the half-filled sys- 
tem. Applying gauge transformations to right-moving 



and left-moving chiral fermion fields independently, we 
can form another duality transformation 

IpRXa = e^'^RXa, ^LXa = ^"^LA., (55) 

which gives the following relations between order param- 
eters 

OcdwM — -O p dw(tt), OdDwM ~ O/dwWi 

Osdw(tt) ~ — O p sdw(tt), OdSDw(Tr) ~ 0/sdw(tt), 

(56) 

where we have omitted fci?A-dependcnt prcfactors. We 
thus find that the transformation (55) exchanges rung- 
centered orders (s- and d-wave orders) and plaquette- 
centered orders (p- and /-wave orders). See Fig. 1. The 
low-energy effective Hamiltonian can be decoupled into 
self-dual parts and symmetry breaking perturbations, in 
which the duality relation is visible. The duality relation 
III thereby gives a parameter mapping as similar to the 
duality relations I and II. In terms of the bosonic fields, 
this transformation is given by the translation of the 4> c + 
field by ir/2 and one can see that the duality relation 
(56) was transparent in the results of the bosonization 
studies. 12,13 We note that this transformation can be 
written in such a compact form only in the low-energy 
effective theory 

6.3 Ising duality 

Finally we mention the Ising duality transformation, 10 
which interchanges the S _ and # s _ fields, 

9 s -^(p s -, (p s -^6 s -, (57) 

and leaves the rest of the <p and 9 fields unchanged. This 
transformation converts a cZDW state into a d-Mott state 
and a CDW state into an s-Mott state. Thus, this is a 
duality between ordered and disordered phases and qual- 
itatively different from the duality relations I-III. 

7. Strong coupling approach at half-filling 

In this section, we discuss the stability of ordered 
states at half-filling using strong-coupling expansions as 
was discussed in ref. 13. Using the duality transforma- 
tions I and II, we can easily find various transition lines 
between ordered (density-wave) and disordered (Mott- 
insulating) phases. 

Let us discuss CDW, dDW, SDW, and dSDW orders. 
For individual phases, an effective Hamiltonian can be 
derived on doubly degenerate rung basis |+) and |— ) in 
the form 

H cff = ]T(^t/t; +1 + hif), (58) 

3 

where r™ (a — x, y, z) denotes the Pauli matrices, and 
Tj |±) = ±|±). This Hamiltonian suggests that doubly 
degenerate ordered ground states appear in K > \h\ and 
disordered states in K < \h\. These transitions are char- 
acterized by the Ising criticality (c = 1/2). 
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In the following derivation, we use convenient basis, 

|i> H,i,t c }aJ >' l 2 > H,u c k T l°>' . . 

(59) 

l 3 >= c 5,i,T c 5,i,J >' l 4 )^ c kT c l2,il°)- 

7.1 CDW-s-Mott 

The perturbation expansion for a CDW ordered state 
can be performed in the strong coupling limit 2Vj_ — 
U — > oo and — U — J±/2 — 71/4 + Vj_ — > oo, where two 
degenerate ground-state basis per rung are |+) = |3) and 
|— ) = |4). By treating hopping terms (tn,t±,t pa i T ) as the 
perturbation, the couplings of the effective Hamiltonian 
are derived as 

K = 2t\f/(2V±-U), (60a) 
h = t pair - 2t ± 2 /(-U + V±- J±/2 - 71/4). (60b) 

A CDW ordered state appears for K > \h\ and an s-Mott 
state for K < \h\. 

7.2 dDW-d-Mott 

Using the duality transformations, we can derive the 
effective Hamiltonian for other ordered states. For the 
dUW phase, the unperturbed ground-state basis are 
|±) = ±{(|1) - |2» T »(|3) - |4))} in the strong cou- 
pling limit V± + tpair + 7i_/2 + 71/4 -> oo and U-V± + 
3t P air + 7j_/2 + 71/4 — > oo. The couplings of the effective 
Hamiltonian are given by 

K = 8t||7(4V ± + 4t pair + 27_l + 71), (61a) 

/l = (-C/ + V± + tpair)/2 - (2J ± + JD/8 

- At ± 2 /(U -V± + 3t pair + J ± /2 + 71/4). (61b) 

One can say that a dDW ordered state appears for K > 
|/i| and a rf-Mott state for K < \h\. 

7.3 SDW-d-Mott 

For the SDW phase, |+) = |1) and |-) = -|2) in the 
strong coupling limit 2U + J]_ — > oo and [/ — Vj_ + t pa i r + 
71/4 — » oo. The couplings of the effective Hamiltonian 
are derived as 

if = 4<||7(7i + 2[/), (62a) 

ft = - J±/2 - 2i ± 2 /(t / ~V± + Vir + 71/4). (62b) 

One can say that an SDW ordered state appears for K > 
\h\ and a rf-Mott state for K < \h\. 

74 dSDW-s-Mott 

For the dSDW phase, |±) = |{(|3) + |4»±i(|l) + |2))} 
in the strong coupling limit V± — i pa ; r — J± /2 + 71/4 — > oo 
and — J7 + Vx —t pa i r — 37^/2 — 71/4 — ► oo. The couplings 
of the effective Hamiltonian are derived as 

K = 8t|| 2 /(4Vx - 4t pair - 27x + 71), (63a) 
h = (U - V ± + £ pair )/2 - (27x - 7D/8 

- 4t ± 2 /(-U + V ± - Vir - 37x/2 - 71/4). (63b) 

One can say that a dSDW (spin current) ordered state 
appears for K > \h\ and an s-Mott state for K < \h\. 
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8. Discussions 

In summary, we have established duality relations 
in correlated electron systems on the two-leg ladder. 
Our arguments clarify mutual relations between conven- 
tional and various unconventional density- wave orders, 
as shown in Fig. 1. The present duality arguments re- 
veal that the stability of unconventional density-wave 
orders such as staggered flux (7DW) and circulating 
spin current (dSDW) is equal to that of conventional 
density- wave orders in the dual parameter spaces. Re- 
cently, large-scale numerical analyses 14 reported the ap- 
pearance of incommensurate 7DW (quasi-)long range or- 
der upon doping. Applying the duality relations to this 
result, we can immediately conclude that new incom- 
mensurate CDW, SDW, and dSDW phases stably exist 
in the doped generalized Hubbard ladder in the dual pa- 
rameter spaces. These duality arguments can be easily 
generalized to the Hubbard ladder including inter-rung 
interactions. 

Next, we found a duality structure in the SO(5) sym- 
metric Hubbard ladder system. The SO(5) symmetry, 
which was proposed to unify AFM and ciSC, also unifies 
(iSDW and sSC. This gives a new route to s-wave su- 
perconductivity in strongly correlated electron systems. 
If the coupling parameters of a half-filled system is close 
to the SO (5) symmetric region and the circulating spin 
current correlation is dominant in the ground state, hole 
doping can be a symmetry breaking perturbation which 
enhances the s-wave superconductivity. Thus this system 
would show a crossover from a g?SDW dominant state to 
a sSC one upon doping. 

Finally, we give a remark on the spin-chirality dual- 
ity which was previously introduced for the spin lad- 
der. 28,29 In the spin-chirality duality transformation, an- 
tifcrromagnctic spin and vector chirality degrees of free- 
dom are converted to each other. One can see that this 
spin-chirality duality has an analogy with the duality I 
for electron systems: in the duality I, SDW (antiferro- 
magnetic spin order) is related to spin current, which is 
expected to have a spin vector chiral order. Actually, Le- 
cheminant and Totsuka 30 showed in a Majorana fermion 
representation that the spin-chirality duality transforma- 
tion can be written as a gauge transformation for the 
fermions, which is similar to the transformation given in 
the present paper. 
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